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We report an improved measurement of the B + — > D + D° and B° — > D°D° decays based on 
657 x 10 6 BB events collected with the Belle detector at KEKB. We measure the branching fraction 
and charge asymmetry for the B + -> D + D° decay: B(B + -> D + D°) = (3.85 ± 0.31 ± 0.38) x 10~ 4 
and A C p{B + -> L>+L>°) = 0.00 ± 0.08 ± 0.02, where the first error is statistical and the second is 
systematic. We also set the upper limit for the B° -> D°D° decay: B(B° -» D°D°) < 0.43 x 10~ 4 
at 90% CL. 

PACS numbers: 13.25.Hw, 14.40.Lb 



Recently, evidence of direct CP violation in B° — > 
D + D~ decays was observed by the Belle collabora- 
tion [l[ , while BaBar measured an asymmetry consistent 
with zero 0- A large direct CP asymmetry for this de- 
cay can indicate the presence of new-physics effects in the 
electroweak penguin sector Q . A similar effect might be 
observable in the charged mode B + — ► D + D°. This de- 
cay has already been observed by Belle 4] and confirmed 
by BaBar @. The latter analysis also includes charge 
asymmetry measurement. 

In this paper, we report an improved measurement of 
the branching fraction and charge asymmetry for B + — > 
D + D° decay and a search for the decay B° — ► D°D°. 
The latter can only be produced by a W exchange di- 
agram. We use a data sample of (657 ± 9) x 10 6 BB 
events collected with the Belle detector at the KEKB 
collider [§]. The inclusion of charge-conjugate states is 
implicit throughout this paper. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cherenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) com- 
prised of CsI(Tl) crystals located inside a superconduct- 
ing solenoid coil that provides a 1.5 T magnetic field. An 
iron flux-return located outside the coil is instrumented 
to detect Kl mesons and to identify muons (KLM) . The 
detector is described in detail elsewhere [7|. For the first 
sample of 152 million BB pairs, a 2.0 cm radius beam 
pipe and a 3-layer silicon vertex detector were used; for 
the latter 505 million BB pairs, a 1.5 cm radius beam 
pipe, a 4-layer silicon detector and a small-cell inner drift 
chamber were used [8j. 

Each track's transverse momentum with respect to the 
beam axis is required to be greater than 0.075 GeV/c 
in order to reduce the combinatorial background. For 
charged particle identification (PID), the measurement 
of the specific ionization (dE/dx) in the CDC, and sig- 
nals from the TOF and the ACC are used. Charged kaons 
are selected with PID criteria that have an efficiency of 
88% with a pion misidentification probability of 8%. All 
charged tracks that are consistent with a pion hypothe- 
sis and that are not positively identified as electrons are 
treated as pion candidates. 



Neutral kaons are reconstructed in the decay Ks — > 
7T + 7r~; no PID requirements are applied for the daugh- 
ter pions. The two-pion invariant mass is required to be 
within 9 MeV/c 2 (~ 3a) of the K° mass and the dis- 
placement of the tt + tt~ vertex from the interaction point 
(IP) in the transverse (r — ip) plane is required to be be- 
tween 0.2 cm and 20 cm. The Ks momentum and the 
vector from the IP to the 7r + 7r _ vertex are required to be 
collinear in the r — ip plane to within 0.2 radians. 

Photon candidates are selected from ECL showers not 
associated with charged tracks. An energy deposition of 
at least 75 MeV and a photon-like shape of the shower 
are required for each candidate. A pair of photons with 
an invariant mass within 12 MeV/c 2 (<~ 2.5<j) of the ir° 
mass is considered as a ir° candidate. We require that 
the 7T° momentum be greater than 0.35 GeV/c in order 
to reduce the combinatorial background. 

We reconstruct D° mesons in the K + tt~ , K + ir~ir + ir~ 
and K + ir~ir° decay channels. The D + candidates are 
reconstructed in the K~ir + ir + and Kstt + final states. 
We require the invariant mass of the D° and D + can- 
didates to be within 11 MeV/c 2 (1.5cr for K + ir~ir° and 
2.5<7 for other modes) of their nominal mass. We perform 
a mass-constrained fit for D candidates to improve their 
momentum resolution. 

To suppress the large background from B + — > D+D° 
with the K + from the Df decay misidentified as a pion, 
none of the pions from D + should be consistent with the 
kaon hypothesis. This requirement has an efficiency of 
93% and kaon misidentification probability of 9%. 

We combine D° and D + (D°) candidates to 
form B + (B°) candidates. These are identi- 
fied by their center-of-mass (CM) energy difference, 
AE = (JV Ei) — -Ebcam, and the beam constrained mass, 
Mbc = v^bcam - (EiPi) 2 > where E hcam is the beam en- 
ergy and pi and Ei are the momenta and energies of the 
decay products of the B meson in the CM frame. We se- 
lect events with M bc > 5.2 GeV/c 2 and \AE\ < 0.3 GeV, 
and define a B signal region of \AE\ < 0.02 GeV, 
5.273 GeV/c 2 < M hc < 5.287 GeV/c 2 . In an event 
with more than one B candidate, we choose the one with 
smallest x 2 from the D mass-constrained fit. We use 
Monte Carlo (MC) simulation to model the response of 
the detector and determine the efficiency @. 

Variables that characterize the event topology are used 
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to suppress background from the jet-like e + e~ — > qq con- 
tinuum process. We require | cos# t hr| < 0-8, where # t hr is 
the angle between the thrust axis of the B candidate and 
that of the rest of the event; this condition rejects 77% 
of the continuum background while retaining 78% of the 
signal. To suppress high background in the D°D° final 
state, we use a Fisher discriminant, T, that is based on 
the production angle of the B candidate, the angle of the 
B candidate thrust axis with respect to the beam axis, 
and nine parameters that characterize the momentum 
flow in the event relative to the B candidate thrust axis 
in the CM frame We impose a requirement on T 

that rejects 52% of the remaining continuum background 
and retains 86% of the signal. 

We determine the signal yield from the two- 
dimensional unbinned extended maximum likelihood fit 
(2D) to the AE-Mbc distribution. The signal proba- 
bility density function (PDF) is described by a dou- 
ble Gaussian for AE and a single Gaussian for Mb C) 
taking into account AE-Mh c correlations. We use the 
B + — > DfD° events in our data sample to calibrate the 
means and resolutions of the signal shape. The contin- 
uum, BB and B + — > DfD° background contributions 
are described separately. We use a linear function for 
AE and a threshold function for Mbc [H| to describe 
the continuum PDF. The BB background is modeled 
by a quadratic polynomial for AE, a threshold function 
for Mbc combined with a small peaking component (a 
wide Gaussian for AE and a Gaussian for M^ c ). The 
shape of the peaking background and parameters of the 
threshold function are fixed from the generic BB MC. 
The AE linear slope and quadratic term are free param- 
eters. The peak in the AE distribution near —70 MeV 
coming from the B + — > D+D° decay is described by a 
Gaussian for AE and a Gaussian for M^ c . Again, we use 
B + -» D+D° to obtain the parameters of this PDF. The 
region AE < —0.1 GeV is excluded from the fit to avoid 
contributions from B — > DD* decays. 

As a cross-check, we also perform separate one- 
dimensional fits to the AE and M^ c distributions, re- 
quiring the other variable to be in the signal region. The 
results are given in Table HI where the listed efficiencies 
include intermediate branching fractions. The projec- 
tions of the 2D fit result are shown in Figs. [1] and [2 The 
90% confidence level (CL) upper limit for B° D°D° 
signal yield is obtained using the POLE program [l2j 
based on the Feldman-Cousins method [HI . The system- 
atic uncertainty (described later) is taken into account in 
this calculation. 

To calculate the charge asymmetry in the B + — » 
D + D° decay channel, we repeat the fits separately for 
the B+ -> D+D° and B~ -> D~ D° samples. The AE 
distributions for B+ -> D+D° and B~ -> D~D° candi- 
dates are presented in Fig. [3J The fit results are given in 
Table Ufl Using the results of the 2D fits, we obtain the 
charge asymmetry: 
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CP 



N{D-D°) - N(D+D°) 
N(D-D°)+N(D+D°) 



0.00 ±0.08 ±0.02 



We calculate the B + — > D + D° branching fraction sep- 
arately for each D decay channel; the results are consis- 
tent with each other. As an additional check, we apply 
a similar procedure to a decay chain with a similar final 
state: B + — > D+D°. We measure the branching fraction 
B(B+ -> D+D°) = (9.5 ±0.2) x 10~ 3 , where the error is 
statistical only. This is consistent with the world average 
value (10.0 ±1.7) x 10~ 3 The charge asymmetry in 
this final state is consistent with zero: (—0.5 ± 1.5)%. 
We also measure the charge asymmetry for the D + D° 
background events and find a value consistent with zero: 
(-1.4 ±1.3)%. 

Table IIIII shows the sources of the systematic uncer- 
tainty. The errors due to knowledge of D branching 
fractions are taken from Ref. (l4| . The uncertainty in 
the tracking efficiency is estimated using partially recon- 
structed D*+ -> D [Kstt + tt-}tt + decays. The uncer- 
tainty in the PID efficiency is determined from D* + — > 
D°[K~Tr + ]ir + decays. The error in signal yield determi- 
nation is estimated by varying the signal and background 
shapes and fit range. We assume equal production rates 
for B + B~ and B°B° pairs and do not include the uncer- 
tainty related to this assumption in the total systematic 
error. 

The asymmetry measurement contains the following 
systematic errors: tracking efficiency difference for 
(0.013), particle identification efficiency difference for 7r ± 
and K ± (0.004) and signal yield determination (0.015). 
The total systematic uncertainty is 0.02. 

In summary, we report improved measurements of the 
B+ -> D+D° branching fraction B(B+ -> D+D°) = 
(3.85 ± 0.31 ± 0.38) x 10~ 4 . The charge asymmetry 
for this decay is measured to be consistent with zero 
A C p{B+ -> D+D°) = 0.00 ±0.08 ±0.02. We also set an 
upper limit for the B° — > D°D° decay branching frac- 
tion of B(B° -> D°D°) < 0.43 x 10" 4 at 90% CL. These 
results are consistent with our previous results Q and su- 
persede them. Our results are also consistent with BaBar 
measurements [5|. 
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TABLE I: Yields from the AE, M bc and 2D (AE-M bc ) fits, detection efficiencies including intermediate branching fractions, 
and corresponding branching fractions. Upper limits are at the 90% CL. 



Decay channel 


AE yield 


M bc yield 


2D yield 


e, 10~ 4 


B, KT 4 


B° -> D°D° 


-4.5 ±29.7 


5.7 ±28.6 


0.4 ±24.8 (< 46) 


16.4 


< 0.43 


B ± -» D ± D° 


366.4 ±31.8 


376.4 ± 30.7 


369.7 ± 29.4 


14.6 


3.85 ±0.31 ±0.38 




AE, GeV M bc , GeV/c 2 



FIG. 1: AE (a) and M bc (b) distributions for the B + — > D + D candidates. Each distribution is the projection of the signal 
region of the other parameter. Points with errors represent the experimental data, open curves show projections from the 2D 
fits and cross-hatched curves show the BB component only. 




AE, GeV M bc , GeV/c 2 



FIG. 2: AE (a) and M bc (b) distributions for the B° —* D°D° candidates. Each distribution is the projection of the signal 
region of the other parameter. Points with errors represent the experimental data, open curves show projections from the 2D 
fits and cross-hatched curves show the BB component only. 



TABLE II: Charged B meson yields from the AE, M bc and 2D (AE-M bc ) fits. Errors are statistical only. 



Decay channel 


AE yield 


M bc yield 


2D yield 


B+ D+D° 


183.9 ±21.5 


184.4 ±21.4 


184.2 ± 20.4 


B~ -> D~D° 


183.4 ±22.1 


192.5 ± 21.8 


185.4 ±21.0 
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0.05 0.1 0.15 0.2 
AE, GeV 




0.15 0.2 
AE, GeV 



FIG. 3: The AE distribution for (a) B + -> D+ D° and (b) B~ D~ D° . Points with errors represent the experimental data, 
open curves show projections from the 2D fits and cross-hatched curves show the BB component only. 



TABLE III: Sources of systematic uncertainty. 



Source 


B(D°D°) B(D+D°) A C p(D + D°) 


D branching fraction 


7% 


5% 





Tracking 


6% 


6% 


0.013 


PID 


2% 


4% 


0.004 


7T° reconstruction 


2% 


1.3% 





N(BB) 


1.4% 


1.4% 





MC statistics 


1% 


1% 





Signal fit 


4% 


4% 


0.015 


Total 


10.5% 


9.8% 


0.02 
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